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Abstract 

rn"; 

^ . We examine connections between rationality of certain indefinite integrals and equi- 

ps) ! librium of Coulomb charges in the complex plane . 

5h ! 

Oh. 

<^ I. Introduction. Adler-Moser Polynomials 

^ In their 1929 paper J.L.Burchnall and T.W. Chaundy [3] examined the following (apparently 
j> . elementary) question: What condition must be satisfied by two polynomials p{z) and q{z) in one 
OO variable z in order the indefinite integrals 



dz, / ^ dz (1) 



O 

cn 
o 

■ may be rational, i.e. expressible without logarithms. Under the assumption that p and q have no 

o |. common or repeated factors the problem is reduced to finding polynomial solutions of the following 

fH [ differential-recurrence equation 

; p'q - 2p'q + pq = (2) 

S It was shown that solutions to (2) are Adler-Moser polynomials [1]: 

: 

><■ P = Oi, q = e,+u 9o{z) = l, e^{z) = z, deg^, = ^^^^ (3) 

The recurrence relation (2) for the Adler-Moser polynomials (8) can be rewritten as the first-order 
(three term) relation [1], [3] 

Cl^n-l - = (2^ + l)^n (4) 

Considering it as a first order linear differential equation for 9n+i we obtain recursively all solutions 
introducing an integration constant at each step. Thus the nth polynomial depends on n free 
parameters 

On = (^n{z;ti,t2, . . . , tn) 

Polynomiality of all solutions is assured by the rationality of (1). 
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Several first examples of the Adler-Moser polynomials are: 

^0 = 1 

di = Z 

e2 = z^ + 12 (5) 

es^z^ + bt2z^ + hz - ml 

In (5), we set ti = , since ti can be absorbed by the translation z ^ z + ti 

It was observed in [2] that the roots of the consecutive polynomials ^j, ^j+i are the equilibrium 
coordinates of i{i + l)/2 positive and (i + + 2)/2 negative Coulomb charges (with values ±1 
correspondingly) on the plane. Namely, the function 

n(n+l)/2 (n+l)(n+2)/2 n(n+l)/2 (n+l)(n+2)/2 

XI X \a.\yi-yj\- ^ ^ hi\xi-yj\ 

i<j=l i<j=l 1=1 j=l 

has a critical point when Xi, yi are roots of 6m On+i respectively. 

The constructive way to find the Adlcr-Moscr polynomials is to use their determinant repre- 
sentation built by means of the Darboux-Crum transformations of the operator (P/dz^. 

This construction can be demonstrated in the following diagram: 

Hi — LiMi — Mj+iLj+i ifj+i = Lj+iMj+i = Mj+2-^i+2 -f^i+2--- (6) 

where 

^ ^ Oiiz) d 9i-,{z) ^ _ Oi^iiz) d 9i{z) 



0i^i{z)dz 0i{z) ' ' 0i{z) dz0i^i{z) 

and Lq = Mq = d/dz. 

From (6), we derive the intertwining differential operators Di,Ui 

Ui-f^^HiUi, -f-Di = DiHi Di = Mi...M,, Ui = L,...Li (7) 
dz^ dz^ 



This, together with the obvious property DiUi — leads to the Wronskian representation 



connecting Hq = d"^/ dz"^ with if, 
This, together with the obvi 
for the Adler-Moser polynomials. In more details 

0n = COnstPr[?/;i, . . . , ^„], i>'^ = i>n-l: i>l = Z (8) 



where the nth polynomial is (up to a constant factor) the Wronskian W = det | |c?Vi/ '^^1|o<i<n-i,i<j<n- 
The intertwining operators are also expressible in terms of Wronskians. For instance 
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II. Indefinite integrals related to equilibrium of charges 



By analogy with the previous section we pose the following question: when does an equilibrium 
condition for charges of different sign is equivalent to rationality of some indefinite integrals? 

Before addressing this question we write the equilibrium equation for two types of charges in 
terms of polynomials: 

Lemma 1 Let 

J2 QiQj^T^\zi- zj\ (10) 

l<i<j<n+m 

be a real function of n + m complex variables ZiEC,i — l...n + m and 

1, i — 1 . . .n 
—A, i — n + l...m + n 



Qi e M, Qi 
Then E has a critical point at 



Xi, i — 1 . . .n 
yi-n: i ^ n + 1 . . .m + n 

{p.«K:=^.W-2A^^=f^ + AM.);?;^ = (11) 

where p, q are the following polynomials 

n m 

P{z) = Y[{z - Xi), q{z) ^Yl{z-yi) 



i=l i=l 



Proof: Is a calculation. We use the partial-fraction decomposition of -^^^ around singular 
points Xi and yi. 

Now we are in position to find an equivalent condition for p and q in terms of rational integrals. 

Indefinite integrals of the algebraic function do not contain logarithmic terms if the function 
does not have singularities of the type z~^. The following lemma shows when certain functions, 
related to solutions of (11) , are free of simple poles and their integrals are rational. 

Lemma 2 Let p{z) = YYi^ii^^^i)^ — YYiLii^ ^ Vi) '^^ ''^^^ have multiple or common roots. 

Then residues of simple poles of (f^/p^ and p^^^ / q^ vanish: 



qizY^ pizf/^ 
Resz=x — r^TT = 0, i — l...n, ReSz=v — , ,^ = 0, i — l...m 

' p{zy ' , z ^^^y 



iff p and q satisfy (11). 
Both indefinite integrals 

p{zY ' J q{z) 
are rational iff p, q satisfy (11) with A = 1/2, 1, 2 



/ -^l^d^ (12) 
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Proof: Let us factor ^(2;) as p{z) = {z — Xi)P(z). Then the condition ReSz=xi^^^ = imphes 
that 



Since P{xi) —p'{xi), P'{xi) —p"{xi)/2 we get the following equation 

p"{xi)q{xi) - 2Ap'{xi)q{xi) = 

Thus, since Xi are roots of ^(2;) 

p"{z)q{z) - 2Ap'{z)q'{z) = a{z)p{z) (13) 

where a{z) is a polynomial. 

Similarly, we get from the condition ReSz=yi ^g^^)2 = that 

AM^)q"{z) - 2Kp\z)q'{z) = h{z)q{z) (14) 

Adding h?'p{z)q" [z) to (13) and p"q{z) to (14) we get 

p"{z)q{z) - 2Ap'(z)q'(z)+AMz)q"(z) = A(z)p{z) = B(z)q{z) 

where A{z) = a{z) + A2g"(^), 5(2;) = b{z) +p"{z). 

Since ^(2;) and q{z) do not have common roots it follows from 

A{z)p{z) = B{z)q{z) 

that A(z) = N{z)q{z), B(z) — N{z)p{z), where N{z) is a polynomial. Thus, we arrive at 

p"{z)q{z) - 2Ap'{z)q'{z) + A'^p{z)q"{z) = N{z)p{z)q{z) (15) 

The degree of the polynomial in the left hand side of (15) is at most deg(g') +deg(p) — 2 while the 
degree of the right hand side is at least dcg(g) + deg(p), unless N(z) = 0. Therefore, N{z) = and 
we obtain (11). The choice A = 1/2, A = 1, A = 2 is essential for the rationality of (12) because 
the both exponents 2 A and 2/A are integer. This completes the proof. 

Since the cases A = 1/2 and A = 2 are equivalent (one is connected with the other by 
permutation of p with q) we have two second order bilinear differential relations corresponding to 
two types of rational integrals 



• Adler-Moser polynomials, A = 1 (see 

{P, ?}i = 'Wqiz) - 2^^ + p(.)^ = 



/ p(2)2 dz, J ^(^)2 dz 



dz^ ' dz dz dz^ 

(16) 
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• A = 2, 



{P, ?}2 = '^q{z) - + 4p(.)^ = 

(17) 

Equation {p, q}2 = has been mentioned in [2] in connection with the problem of equihbrium 
of the point vortices in two-dimensional hydrodynamics. 

We devote the sequel of the work to investigation of the second case and its generalizations. 



III. A = 2 

We begin by summing up what we have obtained so far. 

Let p{z) and q{z) be polynomials of the complex variable z which do not have common or 
multiple roots. 

We study the following three problems: 



When the indefinite integrals 



are rational ? 



-^dz, I — ( 



dz (18) 



• When the system of n and m Coulomb charges of values 1 and —2 respectively can be in 
equilibrium in the complex plane and what are coordinates of the charges ? 

In other words when does the energy function 

E = J2Q,Q,log\z,-z,\, Q.= { i:2/7=n + l..n + m ^^^^ 

have a critical point and at which Zi,i — l..n + m ? 

• When does the equation 

{p, q}, = (20) 

have polynomial solutions? 

The answer to questions (18), (19) and (20) is given by 

Proposition 1 Let i & 'Z, and (modulo translation of z) p_i = z, Po = qo = 1, Qi — z. Then 
all solutions of the recurrence relations 



p'iPi-i - PiP'i-i = (6^ - ^)qf 



(21) 
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are polynomials of the degrees 

deg{pi) = i(3i + 2), deg(gi) = i(3i - l)/2, i e Z 
T/ie?/ provide all solutions of (20) 

{Ph (li\2 = {Ph ?i+i}2 = 0> ^ e Z 
Integrals (18) are rational iff p and q are connected by the bilinear equation (20) i.e. iff 

P^Pi,(l^ Qi, or p^pi,q^ 

The energy function (19) has a critical point provided Zi are zeros of above pairs of polynomials: 
zeros of p, and q being positions of charges 1 and -2 respectively. 

Proof: 

• The equivalence between the bihnear equation (20) and rationahty of integrals (18) is a 
corollary of Lemma 2. 

• The equivalence between (20) and existence of critical points of the energy function (19) is 
a corollary of Lemma 1. 

• Since p — z^ -\- . . . , q = z'^ + . . from the highest symbol of (20) we get the Diophantine 
equation connecting n with m 

(n — 2m)^ — n + Am = 

Its solutions are 

n — ni,m — mi, or n — ni,m — rrii+i, ieZ 

where 

ni^i{3i + 2), rui^ ^ 

Let {pi, qi}2 = 0. Considering this equation as a second order differential equation with a 
solution 5j, by elementary methods we find that its second linearly independent solution is 
given by 

/Pi 
-^dz, {pi, qi+i}2 = (22) 

and is polynomial by the rationality of (18). The degrees of polynomials are connected by 

deg(gi+i) + deg(gi) = deg(pi) + 1 

It is seen that this relation is satisfied if deg(pi) = Ui, deg(gi) = m^. 
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The similar procedure holds if we fix qi and consider and pi as linearly independent 
solutions of (20): 



By freedom in choosing linearly independent solutions of (20) we write analogs of (22) and 
(23) for decreasing indices 



and we can generate pi, qi by induction in either directions starting at some i. 

Since and are strictly increasing for i > (strictly decreasing for i < 0) and mo = 
no = 0, the induction terminates for these two branches dX p^ = q^ = 1. 

Rewriting (22), (23) or (24) in the differential form, we get the first-order recurrence relations 
(21), which completes the proof. 

Here we write down several examples of polynomials satisfying above conditions for the branch 




(23) 




(24) 



i > 



go = 1 Po 
qi^ z pi 
q2^ z^ + T2Z - 4ii p2 




and for the branch i < 



Po = 1 
P-i = z 



go = 1 

q-i ^z^ + T_i 



P-2 = + f T-i-z^ + Url^z'^ 
+28x3^^2 + t-2Z - 7t1i 



q_2 = 2;^ + 7t^iZ^ + 35tIiZ^ 

+T_iT_2 - 1^-2 



+ T^2Z- 



:2 - 35t!iZ 



where U and Tj stand for arbitrary constants. 



IV. Intertwining and factorization 



One can try to find polynomials for A = 2 explicitly by analogy with the Adler-Moser case. 
It easy to observe that equation (20) can be written in two different Scrodinger forms 



connected by the following factorization 



' ' dz q ^ dz p 

Since p = and p = Pi are solutions of (20) with q = qi, we can put any of them in (25 a) with 
q = qi- The similar fact holds for (25 b), with q ~ q^ or q = gj+i and p = Pi. This ambiguity in 
the choice of solutions results in a possibihty of different factorizations for the same second order 
operator. We remind that a similar possibility led to explicit representation of the Adler-Moser 
polynomials for A = 1. 

Now, the factorization diagram has the following form 



i 



1i 



(26) 



where 



dz qi qf dzpi-i dz qi qf dzpi 

The horizontal arrows in (26) correspond to factorizations connected with different choices of p at 
fixed q and vice versa, while the vertical arrows correspond to permutation of factors. Due to the 
change of (prefactors of) the second order operators in the horizontal lines we do not have simple 
intertwining relations between 5*0 = jdz^ and Si. 

Note that now there arc two different U^,U^ (left, right) intertwining operators instead of one. 
This makes the line of approach used for Adler-Moser polynomials unpromising for A = 2. 

V. Baker- Akhieser functions, charges in homogeneous field 

Again we ask the question: When the following indefinite integrals 

' eyiY)ikz)dz, exp(2A;^) / — — — e-xjpi—2kz)dz (27) 
q{k,zY J p[k,zy 

are rational? In (27), k is any complex number. 
We call the function 

*(^.-) = ^^exp(fc.) (28) 

rational Backer- Akhieser function for A = 2 by analogy with A = 1 (see below). 
Similarly to the section III we can formulate the following 
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Proposition 2 Let polynomials p{k, z),q{k, z) in z do not have repeated or common roots. The 
following three statements are equivalent: 

• The indefinite integrals 

exp{-kz) J '^{k,z)dz, exp(2A;2;) j -^{k^zy^dz 

are rational in z functions, expressible without logarithms 

• The energy function 

E = kY,i QiZi + Y.i<j QiQj In l^;^ - Zj\, 



= i 'o ^ • ^ , > n = deg(p),m = deg(g) 

[ —2, i — n + l...m + n ovr-/. ov^/ 

has a critical point at 

Xi{k), i — 1 . . .n 



yi-n{k), i ^n + l...m + n 
where Xi,yi are roots of 

n m 

P^Yl{z-Xi{k)), q^Yl{z -yi{k)) 

i=l 

• p,q satisfy the following equation 



1=1 i=i 



(29) 



= {p, q}r ■■= {P, q}2 + 2k{p'q - 2q'p) (30) 



Proof: Follows arguments similar to Lemmas 1, 2. 

One can easily see that equation (29) describes equilibrium of two types of charges with values 
1, -2 in the homogeneous electric field. 

Another observation consists in the fact that 

p{k,z)=p{C), q{k,z) = q{C), C = kz (31) 

Indeed, one can easily check this fact by scaling Zi,i — l...n + m in (29), or by change the 
integration variable z ^ ( = kz in (27). 
Thus, we can set /c = 1 whenever k ^ 0. 

Equating highest symbols in (30) we obtain the following simple 
Lemma 3 Let p, q satisfy conditions of proposition 2 then their degrees are related by 

deg{p) — n — 2m — 2 deg(g) 
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In other words , the total charge has to be zero in order the system to be at rest in the 
homogeneous field. 

Although the case A = 2 has some similarity with A = 1, the same approach to explicit 
representation of Pi, qi , as shown in the section IV, does not apply any more. 
Here we write down several first examples of Pi{k, z), qi{k, z) 

qo{k,z)^l, po{k,z)^l, 

q,ik,z)=C: Pi{Kz) = C-3C + ^: 

q,{k,z)^e+he+'-^c, 

P2{k, z) = C+ (-9 + 2 12) C + (40 -15t2 + 5/3 ts') C' 

+ (-96 + 52t2- 10^2' + 2/3^2') + (112 - 90 12 + f ^2' - 3t2' + 1/9^2^) C' 



(32) 



+ (-48 + 66 - 28 + 5 ^2^ - 1/3 ^2^) C - 18 ^2 + 10 ^2' - 3 ^2^ + 1/3 ^2^ + 48 

where t2 is an arbitrary constant. 

It is easy to show that solution (32) is generic for m = 3 , in the sense that it is not contained 
in a wider class. In other words , coefficients of (properly normalized) q2 cannot depend on more 
than one free parameter. Indeed , we can always take q2 to be monic and without constant term 
(by fixing one root at ^ = 0). If q2 depended on more than one parameter, then it would be 
q2 = z^ + t2Z^ + Tz, where t2 and T are arbitrary constants. In such a case (30) become a linear 
equations for the coefficients of P2- It is verified by elementary linear algebra that this system is 
incompatible. Thus q2 can depend only on one free parameter (say 12). 

It is difficult to get examples of polynomials of higher degrees, even with the help of the 
computer. It is also nontrivial problem to define degrees of q{k^ z). If one assumes that by analogy 
with A = 1 case the order of qi{k, z) is a quadratic function of i (which certainly holds for = 0, 
see proposition 1), then it follows from the above examples that deg(gi(A;, z)) = i{i + l)/2. 

VI. Scaling, reparametrization and bispectrality 

A statement similar to proposition 2 holds for the case A = 1 with the Baker- Akhieser function 
■0 = {p/q) exp(A;2;) instead of (28) and symmetric integration conditions 



The electrostatic analogy describes charges 1, -1 in the homogeneous field and equation (30) is 
replaced by 

= {p, q]i + 2k{p'q - q'p) 
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This equation can be rewritten in the Schrodinger form 

mk,z)=(^-^ + V^i;{k,z)^k'^{k,z), y = 2^1ogg (33) 

It follows from (7), (9) (for more details see [1]) that the solution of (33) is obtained by action of 
the intertwining operator (9) to the eigenfunction ex'p{kz) of the free Schrodinger operator (f/dz^. 
This solution is the Backer- Akhieser function: 

tjj = W[ijji, ...,tpn, exp{kz)]/W[ijji, ...,tpn] 

Therefore, according to (8) 

p{k,z) ^ ex.p{-kz)W[ilJi, . . . ,ijjn,ewikz)], q{k,z) = . . . , V'n] = const^„(2;) 

On the other hand , by arguments leading to (31) , we see that 

^,+V^m = m, V^2-^\ogq (34) 
with 

m -^,^M0, P(C) =exp(-C)iy[V'i(C),---,V'n(C),exp(C)], g(C) =const^„(C) 

We observe that potential V = V{z) in (33) is independent of A;, while potential V = V^(C)- in 
(34) , depends on k. This is due to the fact that rescaling of variable z ^ ( = kz can be absorbed 
by appropriate changes of the parameters (see e.g. (5)) 

On{z; h, t2, ...,tr,) = k-^^+'^^'Hnikz- kt,, kH^, A;^-^^) 

since the Adler-Moser polynomials can be seen as homogeneous polynomials of variables z and 
> with weights 1 and 2i — 1 correspondingly. 

Thus the problem of equilibrium of charges 1 and —1 can be reduced to the spectral problem 
for the fc-independent operator L (33). 

Similarly, one can ask if ^(/c, z) is simultaneously a solution of a spectral problem for z inde- 
pendent differential with respect to k operator A 

AiP{k, z) = Q{z)ilj{k, z), dA/dz = 

where Q{z) is a function of z. 

It turns out [5] that such operators exist and belong to the so called odd bispectral family. 

In general, the problem of finding functions satisfying simultaneously a differential equation in 
z with /c-dependent eigenvalues and a differential equation in k with 2;-dependent eigenvalues is 
called the bispectral problem [5] . 



11 



Now returning to the main subject of this work , one might expect that the A = 2 case can 
belong to some other bispectral (e.g. even [5]) family of the differential operators. Unfortunately , 
the A = 2 case is not , in general , related to the bispectral problem for a second order differential 
operator. Indeed , equation (30) can be rewritten in the following two forms 

(^ + «(C)) ^(0 = ^(0, u{C) = G^logq{0 

(^+^(c))^(c)-^/' = ^(c)-'/^ ^(c) = !^iogp(c) 

We use (32) as a counterexample to the bispectrahty. It is seen that there is no change of 
parameter t2 — > f{ti,t2) combined with the translation z ^ z + ti in (32) , such that p2 or 
q2 become homogeneous (with some weights) in z and t2- Thus, in difference from the A = 1 
case, rescaling of variable z cannot be absorbed by a change of parameter t2. As a consequence , 
potentials u, v in (35) cannot be 2; or A; independent. Therefore , it is impossible , in general , to 
set a bispectral problem with the second order differential operator in z (or k) when A = 2. 



VII. Integrable dynamics of charges 

The charge configurations studied above can be considered as fixed points of some dynamical 
systems [7]. 

Lemma 4 Let two polynomials p{z,t) — YYi=i{^ ~ -^j(^)) ^'^^ Qiz^t) — YYi=^i{^ ~ -^i(O) satisfy 
the following bilinear equation 

Then the roots Zi{t),i = 1 . . .n + m satisfy the system of ordinary differential equations 

dt t Qi dzi 

where 

^ r 1, i^l...n 

\ -A, i^n + l...m + n 

It follows from the above lemma that the critical points of the energy (10) are fixed points of 
(37). Although, as seen from the previous considerations, existence of such critical points depends 
on values m, n and A, equation (36) has polynomial solutions for any integer m,n >0 and real A. 

Although (37) (or equivalently (36) ) is not a Hamiltonian system, it can be embedded in the 
Hamiltonian fiow: 
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Proposition 3 (37) is a trajectory of a system with the Hamiltonian 



Tj ( dzi\'^ ^ QiQjiQi + Qj) /OS 

1=1 ^ ' t<J=l ^ ■" 



In other words, the Hamiltonian equations of motion 

J=f=t=l ■' 

are corollaries of (37) 

Proof: This proposition is a corollary of more general result proved in [7] . 

This is a nontrivial fact that the Hamiltonian system of the point masses interacting through 
two body forces is a corollary of equations of motion of a lower order for pairwise interacting 
points. In general , such a property is equivalent to compatibility of a highly overdetermined 
system of equations. Perhaps , elliptic generalizations of (37) are the only systems satisfying such 
conditions [7]. 

It is easy to note that, when A = 1 (38) is a sum of two independent Hamiltonians 

1=1 ^ ^ i<j=l ^ i=n+l ^ ^ i<j=n+\ ^ ' ■> ^ 

Each of them belongs to the completely integrable (in the Liouville sense) Calogero-Moser system 
[8]. As a corollary (37) is integrable when A = 1 . 

Equation (38) is also the Calogero Moser Hamiltonian if A = —1. 

In general, (38) is not reduced to known integrable Hamiltonians. 

It is conjectured in [7] that (37) is integrable for any real A in the sense that there exist 
2(n + m) — 1 functionally independent real constants of motion which are rational functions of Zi 
and z* (asterisk denotes the complex conjugation and i — 1 . . .n + m). 

It can be seen from the above the case A = 1 has a particular significance both in the problem 
of the rational integrals (1) and from the point of view of the integrable Hamiltonian systems. It 
is interesting to establish the role of the second case A = 2 in the theory of the dynamical system 
(37). 



VIII Conclusion 

The principal question of this paper to be addressed is about existence and explicit representation 
of the Baker-Akhieser function (or equivalently p(/c, z), g(A;, z)) for A = 2. It is worth to mention a 
useful generalization (interesting by itself) which might help to find the answer: By analogy with 
the A = 1 case, one may consider equilibrium of charges on the cylinder. Remind that (e.g. see 
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[4], [7]), equilibrium configurations for A = 1 are zeros of trigonometric polynomials obtained by 
finite numbers of Darboux transformations from the free Schrodinger operator on circle. In such 
setting, both the Adler-Moser polynomials and A = 1 rational Baker-Akhieser function can be 
obtained choosing different parameters in the limit when the radius of the cylinder goes to infinity. 

It might be interesting to examine a similar procedure getting an analog of proposition 1 on 
the cylinder. 

Another interesting question , as was mentioned before , is to understand the role of the case 
A = 2 in the theory of dynamical systems (37) and integrability of related Hamiltonians. 

Finally , we would like to mention relation with the Reduced Model of Superconductivity 
(Richardson, [9]) , the Gaudin Magnet Model [6] from the one hand and the A = 2 equilibrium 
configurations from the other hand. These models are sets of commuting quantum Hamiltonians 
acting on a finite-dimensional Hilbert space. Equations for the Hamiltonian eigenvalues are equi- 
librium conditions for a set of charges of value —2 located in the complex plane, subject to mutual 
repulsion , and attraction of charges of value 1 located at positions defined by the parameters of 
the Hamiltonian. The A = 2 case considered in the present work corresponds to the Richardson 
(Gaudin) models with special restrictions imposed on the parameters of Hamiltonian.lt is inter- 
esting to examine relations between such models and establish their special properties connected 
with the A = 2 case. 
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